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Dwarf Novae in the Shortest Orbital Period Regime: 
II. WZ Sge Stars as the Missing Population near the Period Minimum 
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Abstract 

WZ Sge-type dwarf novae are characterized by long recurrence times of outbursts (~ 10 yr) and short 
orbital periods 85 min). A significant part of WZ Sge stars may remain undiscovered because of low 
outburst activity. Recently, the observed orbital period distribution of cataclysmic variables (CVs) has 
changed partly because outbursts of new WZ Sge stars have been discovered routinely. Hence, the estima- 
tion of the intrinsic population of WZ Sge stars is important for the study of the population and evolution 
of CVs. In this paper, we present a Bayesian approach to estimate the intrinsic period distribution of dwarf 
novae from observed samples. In this Bayesian model, we assumed a simple relationship between the re- 
currence time and the orbital period which is consistent with observations of WZ Sge stars and other dwarf 
novae. As a result, the minimum orbital period was estimated to be 70 min. The population of WZ Sge 
stars exhibited a spike-like feature at the shortest period regime in the orbital period distribution. These 
features are consistent with the orbital period distribution previously predicted by population synthesis 
studies. We propose that WZ Sge stars and CVs with a low mass-transfer rate are excellent candidates for 
the missing population predicted by the evolution theory of CVs. 
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1. Introduction 

Cataclysmic variables (CVs) are close binary systems 
containing a white dwarf and a Roche-lobe filling normal 
star. They have several subclasses, such as, dwarf no- 
vae (DNe), nova- like variables, novae, and magnetic CVs 
(Warner 1995). CVs are considered to be one of the final 
stages of the evolution of low mass binaries. The forma- 
tion and evolution of low mass binaries can be explored 
through the orbital period (Poib) distribution of CVs (for 
a review, see King 1988). 

Stable mass-transfer processes in CVs are driven by an- 
gular momentum removal from the binaries. In short pe- 
riod CVs having Porb ^ 3 hr, the driving mechanism is 
considered to be angular momentum removal associated 
with gravitational radiation (Paczynski 1981; Rappaport 
et al. 1982). A CV evolves toward a short-Porb region 
by losing angular momentum. This phase is terminated 
when the core of the secondary star becomes degenerate. 
Then, the degeneracy pressure changes the mass-radius 
relationship of the secondary star. As a result, the CV 
evolves toward a long-Porb regime after passing through 
the minimum Porb (Pmin)- 

The Porb distribution near Pmin has received attention, 
since the observed Porb distribution has different char- 
acteristics from that expected from theoretical studies. 
It is well known that Pmin is ~ 80 min in the observed 
CV population (Paczynski 1981; Rappaport et al. 1982). 



However, the theoretically predicted Pmin is 60-70 min, 
significantly shorter than the observed value (e.g. Kolb 
1993). This is the so-called "period minimum problem". 
According to population synthesis studies, most CVs 
have already passed through Pmin- It has been proposed 
that the evolution time-scale is long for an evolutionary- 
advanced system with a low mass secondary star, in other 
words, the mass-transfer rate (M) from the secondary star 
is low (e.g. Kolb 1993; Howeh et al. 1997). As a result, 
a spike-like feature is expected to appear near Pmin in 
the Porb distribution due to the accumulation of systems. 
However, the observed distribution is rather flat. This is 
the so-called "period spike problem" (Kolb, Baraffe 1999; 
Renvoize et al. 2002). 

Thus, it is known that there is a "missing" population 
in the observed CVs with a Porb range of 70-80 min when 
compared with the predicted population. Littlefair et al. 
(2008) have recently reported three CVs having a brown 
dwarf secondary, which are candidates for the missing 
population (Littlefair et al. 2006; Littlefair et al. 2007). 
This discovery indicates that CVs can actually reach the 
post-Pmin regime as predicted by the theories, while the 
number of such systems is still not properly known. 

It has been suggested that the problems can be rec- 
onciled by another mechanism of angular momentum re- 
moval in addition to gravitational radiation (RezzoUa 
et al. 2001), such as the effects of magnetic stellar wind 
braking (King et al. 2002), the circumbinary disk (Willems 
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et al. 2005), and a combination of magnetic propeller and 
accretion disk resonances (Matthews et al. 2006). 

On the other hand, it is possible that the missing pop- 
ulation exists, but remains undiscovered. According to 
the disk instability model for DNe, the recurrence time 
of superout bursts (supercycle; Tg) is longer in a system 
with a lower M (Osaki 1995). WZ Sge-type DNe are 
known to have long T, 10 yr), which indicate a low 
M (HoweU et al. 1995; Kato et al. 2001). It has been 
proposed that several WZ Sge stars, or CVs with a low 
M have remained undiscovered because of their low ac- 
tivity (Meyer- Hofmeister et al. 1998). Since WZ Sge stars 
are found concentrating in a short Porb regime, these ob- 
jects are a candidate for the missing population (Patterson 
1998). However, their intrinsic contribution to the whole 
CV population has been poorly understood because their 
Ts were uncertain and their quiescent luminosity is low. 
Gansicke et al. (2009) have recently reported that the 
CV sample from the Sloan Digital Sky Survey (SDSS) 
database shows a period spike feature in the Porb range of 
80-86 min. The spike feature is a consequence of deep im- 
ages and a homogeneously selected CV sample obtained 
by SDSS, which could detect intrinsically faint CVs with 
a low M. Recently, in addition, Tg of several WZ Sge 
stars have been well determined owing to long monitoring 
of DNe by amateur observers (e.g. Waagen et al. 2007) 

In this paper, we investigate the contribution of long- 
Tg DNe to the whole population of CVs by estimating 
their intrinsic population. In § 2, we demonstrate that the 
whole population of CVs has changed significantly due to 
a recent increase of the number of short Porb CVs. In 
§ 3, we estimate the population of DNe by considering 
the detection probability of superoutbursts using Bayesian 
analysis. We discuss the implications of our results to the 
study of CV evolution in § 5. In the final section, we 
summarize our findings. 

2. Recent change in the observed Porb distribu- 
tion in CVs 

Both the period minimum and spike problems arise 
from the discrepancy between the Porb distribution pre- 
dicted by the population synthesis studies and the ob- 
served one. In previous studies, the predicted distribu- 
tion was occasionally compared with the distribution of 
all known CVs. The left panel of figure 1 shows the Porb 
distribution of all known CVs with Porb < 130 min discov- 
ered before 2003 (black) and 2008 (white). The samples 
were selected from Ritter & Kolb catalog ^: 184 sources 
in RKcat 7.0 for the 2003 sample and 312 sources in RKat 
7.9 for the 2008 sample. Both magnetic and non-magnetic 
CVs are included in the samples. As can be seen in the 
figure, the fraction of short- Porb systems in the 2008 sam- 
ple is apparently larger than that in the 2003 sample. We 
applied the Kolmogorov-Smirnov (KS) test to evaluate 
whether the two samples differ significantly. We calcu- 
lated KS probabilities, Pks, from the 2003 and 2008 sam- 

^ (http://www.mpa-garching.mpg.do/RKcat/} 
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Fig. 1. Poi-b distribution of CVs. Left panel: Sources in the 
2003 sample (RKcatT.O, black histogram) and 2008 sample 
(RKcat7.9, white histogram). Right panel: Normalized cu- 
mulative distribution functions (CDF) of the 2003 (solid line) 
and the new samples discovered between 2003 and 2008 (dot- 
ted line). 



Table 1. Results of Kolmogorov— Smirnov tests between the 
2003 (RKcat ver. 7.0) and 2008 (RKcat ver. 7.9) samples. 



sample 




Nl 


4 


only new objects 


0.07 


184 


128 


all objects 


0.73 


184 


312 



* Pks: Kolmogorov-Smirnov probability. 



t A'^i, N2: the number of objects in the 2003 and 2008 samples. 

pies, and also from the 2003 sample and 128 new sources 
discovered between 2003 and 2008. The results are listed 
in table 1. The test for the 2003 samples and the sample 
of new sources gives a low probability of 0.07, and based 
on this we conclude that the distributions of the samples 
differ significantly. For the test with the full 2008 sam- 
ple, the probability is also low at 0.73 and hence it must 
be concluded that the 2003 and 2008 samples differ over- 
all. The difference between the samples is also clear in 
the normalized cumulative distribution functions (CDF) 
of the 2003 sample and the new sources in the 2008 sam- 
ple, as shown in the right panel of figure 1. These results 
demonstrate that the Porb distribution of all known CVs 
has significantly changed recently. 

This change is mainly caused by the outburst detec- 
tion of new DNe and large spectroscopic surveys. First, 
outbursts of new and known WZ Sge stars have routinely 
been discovered recently and hence the number of known 
WZ Sge stars has rapidly increased (Shaffer et al. 2007; 
Waagen et al. 2007; Templeton et al. 2007; Uemura et al. 
2008; Soejima et al. 2009; Shugarov et al. 2008; Zharikov 
et al. 2008). The dramatic increase of newly recognized 
SU UMa-type and WZ Sge-type DNe is also apparent in 
tables in Kato et al. (2010), in which 14 WZ Sge-type 
DNe (out of a total of 32 systems; Ritter, Kolb 2003; 
Kato et al. 2010) was recognized after 2005. This increase 
of the outburst detection is attributed to the development 
of automatic sky monitor systems (e.g. Pojmanski 2002; 
Burd et al. 2005; and most recently, Drake et al. 2009), 
and improvements in the techniques and instruments of 
amateur observers (notably by K. Itagaki). The devel- 
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opment of online networks for transient phenomena has 
also played an important role in circulating information 
on WZ Sge stars, and enabled timely observation essen- 
tial for identifying their nature (Kato et al. 2004). 

Second, large spectroscopic surveys have also played an 
important role in increasing the number of CVs recently, 
for example, SDSS and the Hamburg Quasar Survey 
(HQS) (Szkody et al. 2002; Bade et al. 1998). Those large 
surveys allow us to investigate the Porb distribution with- 
out complicated observation bias which should be consid- 
ered in a sample containing all CVs discovered by various 
detection methods. The Porb distribution of all known 
CVs would be same as before if the distribution of those 
new CVs is same as that of the previously known CVs. 
However, the distribution of the SDSS CVs has a period 
spike feature, and is different from the previously known 
distributions (Gansicke et al. 2009). 

Thus, it is evident that the period minimum and spike 
problems should be discussed not with the sample contain- 
ing all known CVs, but with the sample in which observa- 
tion bias for the Porb distribution can be simply evaluated. 
In the next section, we collect outburst-selected DN sam- 
ples, and estimate the intrinsic population of DNe without 
the complicated bias for the Porb distribution. 

3. Bayesian estimation of the intrinsic population 
of DNe 

As well as new WZ Sge stars, outbursts of known 
WZ Sge stars have been observed routinely. Thereby, we 
can now estimate Tg more accurately for them. In Uemura 
et al. (2010) (hereafter. Paper 1), we revised Ts of several 
DNe and found that WZ Sge stars form a major group in 
Porb ^ 86 min. This implies that a DN enters the WZ Sge 
regime on an evolutionary course at Porb ^ 86 min as M 
decreases. Using Tg and Porb from Paper 1, we consider a 
statistical experiment to estimate the intrinsic Porb distri- 
bution of DNe from observed systems based on a simple 
model for T^. Using Bayesian analysis, we calculate the 
posterior distributions of parameters that characterize the 
Porb distribution (e.g. Trotta 2008). 

3.1. Overview of analysis 

We first present an outline of the analysis. The models 
and samples are detailed in subsequent sections. 

We consider a DN sample, x = {Porb.i}, whose super- 
outbursts were detected at least once in a certain period 
of time. The intrinsic Porb population of DNe can be re- 
produced from this outburst-selected sample if can be 
described with a function of Porb ■ The probability density 
function (PDF) of observed DNe, (5(Porb), is given by the 
intrinsic period distribution, /(Porb), and the selection ef- 
fects arising from the observation process, £)(Porb), in the 
form Q = D ■ I /Aq, where Aq is a normalization factor. 

With a parameter set, y, that characterizes /, the like- 
lihood of X, can be defined by 



Posterior distributions of parameters, P(y | x), are est! 
mated with Bayes' theorem, as follows: 

i(y|x)P(y) _ Q(x|y)P(y) 



^(y I x) = 



(2) 



Ei(y|x)P(y) EQ(x|y)P(y)' 

In our analysis, we used flat distributions for the prior 
probabilities of parameters: P(y) = const. 

3.2. Sample 

Samples were selected from SU UMa-type DNe in 
RKcat7.9. The "ASAS" sample consists of systems in 
which superoutbursts were recorded at least once in the 
ASAS-3 database between 1 Jan 2003 and 31 Dec 2007 
(Pojmanski 2002). A superoutburst is defined as an 
outburst lasting for > 5 d with a small fading rate (~ 
0.1 magd~^). Porb of our sample was limited to < 130 min 
because our simple model of D cannot reproduce the fea- 
tures of a period gap of 2 hr < Porb ^ 3 hr (for details, 
see the next section). We made another sample, called 
the "mixed" sample, using the procedure for data selec- 
tion as for the ASAS sample, but using both the VSNET 
and ASAS-3 databases for searching superoutbursts (Kato 
et al. 2004). 

We excluded two minor groups of objects in the shortest 
Porb regime. The first group is ER UMa stars (ER UMa, 
IX Dra, RZ LMi, V1159 Ori, and Dl UMa), which are 
characterized by atypically short (20-50 d) despite 
a short Porb (Kato, Kunjaya 1995). The other group 
was identified in Paper 1 for the first time and called 
"Group X" . The members of this group have rather aver- 
age Ts (^ 400 d) despite a short Porb- These two groups 
were excluded from our sample because they may lead to 
an overestimate of the intrinsic population of short-Porb 
systems having long Tg. 

We obtained 42 and 146 sources for the ASAS and 
mixed samples. The ASAS sample is expected to be an 
ideal sample without complicated selection effects since 
the superoutbursts were detected through a single obser- 
vation process defined by the ASAS system. The mixed 
sample may be affected by complicated selection effects. 
We calculated the KS probability of both samples. The 
obtained probability is unexpectedly high, 0.98, which 
means that there is no significant difference between those 
two samples. This implies that our method of data selec- 
tion has the advantage of reducing selection bias. The Porb 
distributions of these two samples are shown in figure 2. 

3.3. Model of the Porb distribution 

The intrinsic Porb distribution, /, must be characterized 
by at least two parameters, namely, one for the degree of 
accumulation near Pmin and the other for the real Pmin- 
We express / by the following formulae: 

-"/P/A/ (p>l) 
{p<l) 



(3) 



P = Porb - Pmin (min) 



(4) 



L(y|x) = Q(x|y)=[]Q(Porb.dy)- 



(1) 



where a is a parameter for the degree of the accumulation 
near Pmin- These are the parameters of our Bayes model, 
that is, y = {a, Pmin}- A large a yields a distribution with 



4 



M. Uemura, et al., 




70 80 90 100 110 120 130 
Orbital period (min) 



0.1 





1 


1 1 1 

• 

• i 


: n=0.2 




SW UMa. • • 














-n=1.0 




• 










GWl-ib 
• 








• 

• 


ri=2.0 


,ri=3.0 


1 1 1 



4 6 8 10 

Orbital period (min) - P^in(=72.6 min) 



Fig. 2. Distribution of Porb of the DN samples. The black 
and white histograms indieate the ASAS and mixed samples, 
respectively. The dashed line indicates the best model PDF 
derived from our analysis (for details, see the text). 

a prominent spike near Pmin- The distribution becomes 
flat in the case of a = 0. For a negative a, the distribution 
has a peak in the long-Porb region. The function / has an 
extremum at Porb ~ ^min + 1- is a normalization factor. 
While this formula of / itself has no physical meaning, it 
can approximate the profile of the Porb distribution below 
the period gap predicted by population synthesis studies, 
with high a (a ~ 1) (e.g. Barker, Kolb 2003). This model 
should be used only below the period gap (2 hr < Porb ^ 
3 hr; Spruit, Ritter 1983), since it cannot reproduce the 
gap structure. 

We considered a simple model for the selection effect, D. 
First, we considered that all the systems are pre-Pmin CVs 
and neglected contributions from post-Pmin CVs. The 
contribution of post-P,nin CVs to our analysis is discussed 
in § 4.2. Second, we used a proposed model for the re- 
lationship between P,. and Porb based on observations. 
In Paper I, we showed that most systems are ordinary 
SU UMa stars in a Porb range of ^ 86 min. The average 
Tg is 470 d in 86 min < Porb < 95 min. On the other hand, 
in the shortest Porb regime of ^ 86 min, the major group 
is WZ Sge stars. Based on these results in Paper I, we 
assumed the following for our model, D: i) Tg is constant 
at 470 d in Porb > -Pcrit- ii) Tg becomes longer with shorter 
Porb in -Porb < -Pcrit- hi) We set Pcrit = 86 min. 

The outburst detection probability, Doutb, is inversely 
proportional to Tg. Taking -Doutb at Tg = 470 d to be unity, 
the above assumptions can be expressed as follows: 



-Doutb (-Porb) 



1 



Pcrit 



(-Porb > Pcrit) 
(-Porb < Pcrit) 



(5) 



The parameter n can be estimated from the observed 
Tg listed in table 5 of Paper I. Using GW Lib (Porb = 
76.7808 min and Tg = 8800 d) and SW UMa (Porb ^ 
81.8136 min and Tg ~ 1200 d), for example, we obtained 
n = 2.5 and Pmin = 72.6 min. We show P^outb calculated 
from the observed Tg of our sample with Pmin = 72.6 min in 
figure 3. The dashed lines in the figure denote the models, 



Fig. 3. Outburst detection probability, P*outbi a- function 
of 

-forb — Pmin- Pmin is assumed to be 72.6 min in the fig- 
ure (for more detail, see the text). The filled circles indicate 
SU UMa stars whose Ts is known. The dashed lines indi- 
cate the models with various n. Two objects, GW Lib and 
SW UMa are labeled. They are used to estimate a lower limit 
of n. 

-Doutb, with various n. This figure shows that the estima- 
tion of the parameters with GW Lib and SW UMa yields 
a lower- limit of n for our sample. As discussed in § 4.1, a 
larger n leads to a larger number of short-Porb DNe. In 
the following analysis, we used n = 2.0, which gives a firm 
lower-limit for n, as can be seen from figure 3. n can be 
smaller than 2.0 only when Pmin > 72.6 min. We check 
whether our estimation of n is consistent with the results 
of the analysis in § 4.1. 

Finally, we considered the relationship between the ab- 
solute magnitude. My, at supermaximum and Porb- Since 
the ASAS sample is a magnitude-limited sample, the se- 
lection bias caused by the Porb dependence of My should 
be corrected. Under an assumption of a uniform space 
distribution of DNe, the number of detectable systems, 
-Dmag! can be written as a function of My as 

A„ag(-Porb)ocl0-°-6^^-(^"-). (6) 

Harrison et al. (2004) reported the relationship between 
My at normal outburst maxima and Porb- In addition to 
their relationship, Pretorius et al. (2007) proposed that 
WZ Sge stars have a constant Afy, independent of Porb- 
We follow their models for My at the supermaxima of or- 
dinary SU UMa stars (My^su) and WZ Sge stars (-/l'/y,wz), 
that is: 



Mv,Su(-Porb) = Ci-0.383Porb 
Mv,Wz(P'orb) =C2 



(-^orb ^ -fcrit) fij\ 



{P. 



orb 



where C'2 is a constant determined by -Afy.su (Pcrit) = 
-A^v.wz (-Pcrit)- The unit of Porb is hour. The constant 
C'l is not important in our analysis. 

Finally, the model of the Porb distribution of DNe is 
expressed by the PDF, Q, using the two parameters as 
follows: 

Q(a,P'min) = Anag ' -Doutb ' /(a, Pmin)MQ (», -Pmin) - (8) 
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Fig. 4. Posterior probability distributions of a and Pmin- 
The solid and dotted lines were calculated from the ASAS 
and mixed samples, respectively. The medians and 68.3 % 
confidence levels of the parameters are also indicated. 

3.^. Bayesian analysis and results 

We calculated the likelihood of the sample using equa- 
tion (1). We produced posterior probability distribu- 
tions of a and P„iin using the Markov Chain Monte Carlo 
method (MCMC) (Metropolis 1953). The calculation pro- 
cedure of MCMC is as follows: At the n-th step of MCMC, 
we obtain the likelihood, L„ = i(a„,Pniin,n | x). We then 
randomly move to another point in the parameter space 
of a and Pmin by adding random values drawn from a 
Gaussian distribution which has zero mean and dispersion 
chosen to efficiently sample the likelihood surface. As a 
result, we obtain Ln+i- We count the parameters of the 
{n + l)-th step if > L„ or L„+i/L„ is larger than a 
uniform random number between and 1. Otherwise, the 
n-th parameters are substituted for the (n-l-l)-th ones. 
After discarding the first 100 steps, we sample every 100 
steps until the size of the sample reaches 10^. This cycle 
forms a set of samples for a and Pmin- We obtained ten 
sets of samples with different initial values of a and Pmin- 
We confirmed that the obtained ten samples had the same 
Porb distribution. We merged them and finally obtained 
the posterior distributions. 

Figure 4 shows the posterior distributions of a (left 
panel) and Pmin (right panel). Each panel contains both 
the results estimated from the ASAS and mixed samples. 
The median and 68.3% confidence level of a and Pmin are 
summarized in table 2 and are also indicated in the fig- 
ure. As expected from the KS test in § 3.2, the parameters 
estimated from the ASAS sample agree with those from 
the mixed sample within the errors. The best model of 
Q for the ASAS sample is indicated by the dashed line 
in figure 2. The intrinsic Porb distribution, /, is shown in 
figure 5. 

As shown in table 2, a was estimated to be high, 
namely, the intrinsic Porb distribution, /, has a spike near 
Pmin, as can be seen in figure 5. Short- Porb systems in 
-Porb < f crit (= 86 min) account for 59 % of the total num- 
ber of samples between Pmin < Porb < 130 min. The es- 
timated Pmin is about 70 min, significantly shorter than 
the observed Pmin of about 78 min. 

The estimated Porb distribution, /, is analogous to those 
expected from population synthesis studies (e.g. Kolb 
1993), in terms of the accumulation of systems near Pmin 
and the short Pmin of ~ 70 min. We suggest that i) DNe 



Table 2. Best parameters derived from Bayesian analysis. 



Sample 


n 


Pcrit 

(min) 


Mv 


a 


P . 
min 

(min) 


ASAS 


2.0 


86.0 


eq.(7) 


0.88+-.- 




all RKcat 


2.0 


86.0 


eq.(7) 




71.I+P7 


ASAS 


1.0 


86.0 


eq.(7) 


o.6i±l!:^^ 




ASAS 


3.0 


86.0 


eq.(7) 


i-07l°:r7 


66.21^.8 


ASAS 


2.0 


82.0 


eq.(7) 


0.601^:1 




ASAS 


2.0 


90.0 


eq.(7) 


1 1 c;-|-0.32 
^■^^-0.36 


69.5l^-^ 


ASAS 


2.0 


86.0 


const. 


64+°- ^3 


68.61^;° 


ASAS 


0.2 


86.0 


const. 


10+°-25 


73. 0^^;^ 



0.1 



0.08 - 



0.06 - 

LL 
Q 
CL 

0.04 - 



0.02 - 




I 1^ 1 1 1 1 1 1 

70 80 90 100 110 120 130 

Orbital period (min) 

Fig. 5. Intrinsic Porb distribution of DNe, I, estimated from 
the observed sample. 

having long Tg are an excellent candidate for the missing 
population near Pmin, and ii) the observed Pmin is sig- 
nificantly longer than the real Pmin because of the low 
outburst activity of DNe near Pmin- 

4. Discussion 

Dependence of a and Pmin on the assumed model 
parameters 

In the previous section, we assumed several parameters 
in our model to estimate a and Pmin, such as n, Pcrit and 
My at supermaximum. Here, we discuss the dependence 
of the result on these parameters. We present a and Pmin 
estimated with several different sets of these parameters 
in table 2. 

As can be seen in table 2, a is sensitive to both n and 
Pcrit- Smaller n or Pcrit leads to a smaller a. a also 
depends on the model of My. It is possible that My 
takes a form different from equation (7). According to 
Harrison et al. (2004), Mv,su for superoutbursts has a 
large dispersion and does not strictly follow the My-Porb 
relationship for normal outbursts. This implies that the 
Porb dependence of My.su is weak. Assuming Mysu = 
Mv.wz — const., we obtained a smaller a as shown in 
table 2. In all the above cases, a takes high values of > 0.6 
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and a spike near Pmin appears in the Porb distribution. 
The spike feature disappears only when n is extremely 
small. For example, we show the case of n = 0.2 in table 2. 
Such a small n is, however, unfavorable for the observed 
^outb--Porb relationship, as shown in figure 3. 

Table 2 indicates that Pmin is sensitive to n; a small 
n yields a large Pmin- On the other hand, the estimated 
^min in table 2 are significantly smaller than the observed 
value, even in the extreme case of n = 0.2. Thus, the 
result is robust within the reasonable ranges of the model 
parameters. 

As mentioned in § 3.3, we used n = 2.0, which is a lower 
limit of n in the case of Pmin = 72.6 min. n can be smaller 
than 2.0 only when Pmin > 72.6 min. As shown in table 2, 
all Pmin satisfies Pmin < 72.6 min, except for the extreme 
case of n = 0.2. Therefore, it is less likely that n is smaller 
than 2.0 in our model. 

Our result highly depends on the assumption for the 
-Doutb~Porb relationship. The assumption might be incor- 
rect if a number of short- Ts systems remain undiscovered 
in the shortest Porb regime, such as Group X identified in 
Paper I. The assumption for the -Doutb^Porb relationship 
should thus be reconsidered if the intrinsic population of 
those objects is found to be much larger than the cur- 
rently observed population. However, it is unlikely that 
there are significant undetected systems because the out- 
burst detection probability must be high in such short-Ts 
systems and hence their outbursts should have already 
been discovered. 

4-.2. Implications for CV evolution 

In this section, we discuss the implications of our result 
to CV evolution, considering the presence of post- Pmin 
CVs, magnetic CVs, and past theoretical and observa- 
tional works in this field. 

In the previous section, we regarded all systems as pre- 
Pmin CVs. However, it has been proposed that about 
70% of CVs have already passed through Pmin (Kolb 1993; 
Howell et al. 1997). It is posible that all known DNe 
are pre-Pmin CVs, while most CVs have already passed 
through Pmin and remained undiscovered. In this case, 
the intrinsic Porb distribution would have a taller spike 
near Pmin than that in figure 5 because of the contribution 
of the undiscovered post- Pmin CVs. 

On the other hand, a fraction of the long-Ts DNe in 
our sample might actually be post- Pmin systems because 
post-Pmin systems presumably have quite low M and long 
Tg. Even if most known WZ Sge stars are post-Pmin CVs, 
the spike structure in figure 5 is still expected to appear 
because the spike feature depends, in principle, on the 
number of long-T^ DNe near Pmin and is independent of 
the nature of the systems, namely, whether they are pre- 
or post-Pmin systems. 

In contrast to the spike feature, the real Pmin is possibly 
sensitive to the number of pre-Pmin systems in the short- 
est Porb regime. If the true fraction of post-Pmin systems 
is large in known WZ Sge stars, then it means that super- 
outbursts of the post-Pmin systems are detectable for us. 
If this is the case, we should be able to detect outbursts of 



DNe below the observed Pmin- The estimated Pmin could 
then be close to the observed value. 

The true intrinsic fraction of post-Pmin systems in the 
shortest Porb regime is poorly known, while recent obser- 
vational studies suggest the presence of post-Pmin systems. 
Townsley & Gansicke (2009) investigated the temperature 
of the white dwarf in CVs as a probe of M. Their result 
shows that there are anomalies having quite low temper- 
atures of the white dwarf compared with most CVs in 
the shortest Porb regime. The M of those anomalies are 
possibly lower than the ordinary ones, and hence, sug- 
gesting that they are post-Pmin systems. Littlefair et al. 
(2008) reported three CVs in which the post-Pmin nature 
was dynamically confirmed (also see, Littlefair et al. 2006; 
Littlefair et al. 2007). Among known DNe, WZ Sge has 
the longest Tg and hence is a good candidate for post-Pmin 
systems. Steeghs et al. (2007), however, reported that its 
secondary mass is so high that WZ Sge is in fact a pre- 
-Pmin system. The nature of WZ Sge implies that most 
DNe whose Tg are known are still in the pre-Pmin regime. 

As mentioned in § 3.3, our model of Doutb assumed a de- 
crease of M toward Pmin in the region Porb < 86 min. This 
was required from the observed T^-Porb relationship, as 
can be seen in figure 3. However, this is apparently incon- 
sistent with theoretical calculations for the CV evolution, 
which predict a rather constant M for pre-Pmin systems 
even near Pmin (e.g- Shatter 1992; Kolb 1993; Barker, 
Kolb 2003). If the theoretical prediction is true, we would 
expect many short-Tg systems in the shortest Porb regime. 
The prediction is, however, unfavorable for the observed 
population, because there are only a few short-Tg systems 
in Porb ^ 86 min, such as Group X in Paper I. This incon- 
sistency implies that the current understanding of the sec- 
ondary star is imperfect, since the theoretical M-Porb re- 
lationship depends on the structure of the secondary star. 
Meyer, Meyer-Hofmeister (1999a) proposed that the mag- 
netic activity of the secondary star reduces rapidly as it 
becomes degenerate, and therby, the viscosity in the disk 
reduces and Tg increases. The observed T^-Porb relation- 
ship may be reprocuded if increases by this mechanism 
in Porb < 86 min. 

It has been reported that the Porb distribution of mag- 
netic CVs is similar to that of non-magnetic CVs (e.g. 
Willems et al. 2005). However, this discussion was based 
on samples obtained from various observation processes, 
such as samples in RKcat. As mentioned in section 2, such 
a sample is not suitable to be compared with the the- 
oretical Porb distribution obtained from population sys- 
thesis studies. Our result suggests that, in the case of 
DNe, the period minimum and spike problems are recon- 
ciled with a number of undiscovered WZ Sge stars. Since 
AM Her stars have no accretion disk, and thereby, expe- 
rience no DN-type outbursts, our result also implies that 
another mechanism to reduce short-Porb systems should 
work for magnetic CVs (e.g. Kolb 1995; Meyer, Meyer- 
Hofmeister 1999b). The reduce of magnetic braking in 
magnetic CVs may lead to a diff'erent Porb distribution 
from non-magnetic CVs (Webbink, Wickramasinghe 2002; 
Araujo-Betancor et al. 2005). 
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Gansicke et al. (2009) have recently reported the discov- 
ery of the period spike feature in the CV sample selected 
from the SDSS database. The appearance of the spike is 
in agreement with our result; we reproduced a spike based 
on observed outburst activity, while Gansicke et al. (2009) 
obtained it directly from the deep survey. They discovered 
that CVs having a white-dwarf dominated spectrum form 
a major population in the CVs making the period spike 
feature. Their spectrum indicates that M of those CVs are 
quite low. Hence, they may be WZ Sge stars, and their 
outbursts could be observed with a long Tg. Actually, 
SDSS J080434.20-K510349.2, a system having a white- 
dwarf dominated spectrum, experienced a WZ Sge-type 
superoutburst in 2005 (Pavlenko et al. 2007). The SDSS 
CVs having a white-dwarf dominated spectrum probably 
have the same nature as the objects that our analysis pre- 
dicts to be present in the shortest Porb regime. 

On the other hand, the implications for the period 
minimum problem are different between our result and 
Gansicke et al. (2009). In our result, the observed Pmin 
is significantly longer than the real Pmin because Ts is so 
long that the detection frequency is low near the real Pmin- 
The SDSS images are so deep that we can expect to detect 
CVs with a low M near the real Pmin- However, a sharp 
cut-off exists at Porb ~ 80 min in the Porb distribution of 
the SDSS sample and only a few sources were found in the 
Porb range of 60-78 min. The period minimum problem, 
hence, remains unresolved in Gansicke et al. (2009). The 
lack of the sources in the Porb range 60-78 min may indi- 
cate that a fraction of known WZ Sge stars are post- Pmin 
CVs, as discussed above. Alternatively, some CVs with a 
very low M may remain undiscovered even in the SDSS 
survey. Gansicke et al. (2009) identified CVs based on 
spectral features of objects that were selected by several 
color criteria. A part of the color criteria plays a role in 
excluding single white dwarfs from the CV sample. The 
criteria, however, possibly excludes evolved CVs which 
have a quite low M, and thereby, have optical spectra 
dominated by the emission from a white dwarf. 

5. Summary 

We estimated the intrinsic Porb distribution of DNe as- 
suming a simple relationship between Tg and Porb- The 
estimated distribution has a spike near Pmin and a shorter 
Pmin than the observed value, which was expected by the 
theory of binary evolution. We suggest that the nature of 
the "missing" population is DNe with long Tg near Pmin- 

We thank B. Gansicke for his valuable comments and 
suggestions on this paper. This work was partly sup- 
ported by a Grand-in- Aid from the Ministry of Education, 
Culture, Sports, Science, and Technology of Japan 
(19740104). 

Appendix 1. Outburst-selected DNe sample 

Table 3 and 4 lists the objects of the ASAS and mixed 
sample used in our analysis described in section 3, re- 
spectively. The object name and Porb were quoated from 



Ritter & Kolb catalog 7.9 (Ritter, Kolb 2003). 
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Table 3. Members of the ASAS sample. 



Object 


Porb (niin) 


Object 


Porb (min) 


Object 


Porb (min) 


Object 


Porb (min) 


GW Lib 


76.8 


J0233-1047 


78.9 


J0137-0912 


79.7 


VI 108 Her 


81.9 


WX Get 


83.9 


J1112-3538 


84.1 


GG Scl 


84.5 


2219+1824 


86.3 


V1040 Gen 


86.8 


AQ Eri 


87.8 


J0918-2942 


88.8 


J1025-1542 


89.3 


V436 Gen 


90.0 


EK TrA 


90.5 


OY Gar 


90.9 


J1600-4846 


91.3 


J1536-0839 


91.6 


J0232-3717 


93.2 


AK Gnc 


93.7 


GO Gom 


94.8 


RZ Sge 


98.3 


DT Oct 


104.5 


V368 Pe 


105.1 


VZ Pyx 


105.6 


GG Gne 


105.9 


J1556-0009 


106.7 


VW Hyi 


107.0 


Z Gha 


107.3 


QW Ser 


107.4 


WX Hyi 


107.7 


RZ Leo 


109.5 


DH Aql 


111.4 


J0219-3045 


112.9 


GU Vel 


113.0 


J0549-4921 


115.5 


V877 Ara 


116.7 


TU Grt 


118.2 


TY PsA 


121.1 


KK Tel 


121.7 


V364 Peg 


122.4 


YZ Gnc 


125.0 


GZ Gnc 


126.9 











Table 4. Members of the mixed sample. 



Object 


Porb (min) 


Object 


Porb (min) 


Object 


Porb (min) 


Object 


Porb (min) 


J0329+1250 


76.0 


GW Lib 


76.8 


J0222+4122 


78.9 


J0233-1047 


78.9 


J0137-0912 


79.7 


SS LMi 


80.2 


J1021+2349 


80.5 


J0025+1217 


80.7 


V455 And 


81.1 


AL Gom 


81.6 


SW UMa 


81.8 


J1839+2604 


81.9 


J1959+2242 


83.5 


WX Get 


83.9 


J1112-3538 


84.1 


GG Scl 


84.5 


FL TrA 


84.8 


J0804+5103 


85.0 


V585 Lyr 


85.4 


2219+1824 


86.3 


J0807+1138 


86.4 


GI UMa 


86.4 


DV Dra 


86.4 


RX Vol 


86.5 


MM Sco 


86.7 


V1040 Gen 


86.8 


V1028 Gyg 


86.9 


UZ Boo 


87.4 


AQ Eri 


87.8 


V1454 Gyg 


87.8 


XZ Eri 


88.1 


V4140 Sgr 


88.5 


J0918-2942 


88.8 


J1025-1542 


89.3 


VI 141 Aql 


89.3 


V402 And 


89.4 


V2051 Oph 


89.9 


V436 Gen 


90.0 


BG UMa 


90.2 


HO Del 


90.3 


EK TrA 


90.5 


TV Grv 


90.6 


J1227+5139 


90.7 


OY Gar 


90.9 


J1600-4846 


91.3 


J1536-0839 


91.6 


MR UMa 


91.6 


Var 79 Peg 


91.7 


DO Vul 


92.2 


J1653+2010 


92.7 


J0232-3717 


93.2 


UV Per 


93.5 


GT Hya 


93.5 


AK Gnc 


93.7 


AO Oct 


94.1 


DM Lyr 


94.2 


GO Gom 


94.8 


V551 Sgr 


94.9 


BG Dor 


95.2 


FQ Mon 


95.8 


SX LMi 


96.8 


V701 Tau 


96.9 


SS UMi 


97.6 


NSV 4838 


97.6 


BZ UMa 


97.9 


J0824+4931 


97.9 


KS UMa 


97.9 


V337 Gyg 


97.9 


RZ Sge 


98.3 


TY Psc 


98.4 


IR Gem 


98.5 


V699 Oph 


98.6 


V1208 Tau 


99.1 


V1504 Gyg 


100.1 


GY UMa 


100.2 


RU Hor 


100.4 


BB Ari 


101.2 


VW GrB 


101.8 


PU Per 


102.7 


FO And 


103.1 


AW Sge 


104.3 


TU Tri 


104.3 


DT Oct 


104.5 


OU Vir 


104.7 


V1212 Tau 


105.1 


V368 Peg 


105.1 


KV And 


105.6 


VZ Pyx 


105.9 


GG Gnc 


105.9 


AY Lyr 


106.1 


lY UMa 


106.4 


J1556-0009 


106.7 


VW Hyi 


107.0 


Z Gha 


107.3 


QW Ser 


107.4 


J0242-2802 


107.4 


BE Oct 


107.6 


WX Hyi 


107.7 


RZ Leo 


109.5 


QY Per 


109.6 


0417+7445 


109.9 


SU UMa 


109.9 


V630 Gyg 


110.0 


J1730+6247 


110.2 


VI 113 Gyg 


110.3 


AW Gem 


110.6 


HS Vir 


110.7 


EG Aqr 


110.9 


TT Boo 


110.9 


V344 Pav 


110.9 


DH Aql 


111.4 


TY Vul 


111.9 


V503 Gyg 


111.9 


PV Per 


112.0 


V359 Gen 


112.2 


V660 Her 


112.7 


BZ Gir 


112.9 


J0219-3045 


112.9 


GU Vel 


113.0 


J2100+0044 


113.8 


NSV 14652 


113.8 


BR Lup 


114.5 


NSV 9923 


114.6 


J2258-0949 


115.2 


J0549-4921 


115.5 


RX Gha 


116.6 


V877 Ara 


116.7 


AB Nor 


117.1 


GP Dra 


117.5 


V369 Peg 


117.6 


TU Grt 


118.2 


HV Aur 


118.5 


EE Peg 


120.5 


TY PsA 


121.1 


BE Ara 


121.7 


KK Tel 


121.7 


V452 Gas 


121.8 


V364 Peg 


122.4 


DV UMa 


123.6 


V419 Lyr 


124.4 


YZ Gnc 


125.0 


DM Dra 


125.3 


V344 Lyr 


126.1 


GZ Gnc 


126.9 


GX Gas 


128.2 


UV Gem 


128.2 
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